Entanglement of formation for an arbitrary two-mode Gaussian state 
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We write the optimal pure-state decomposition of any two-mode Gaussian state and show that its 
entanglement of formation coincides with the Gaussian one. This enables us to develop an insightful 
approach of evaluating the exact entanglement of formation. Its additivity is finally proven. 



PACS numbers: 03.65.Ud, 03.67.Mn, 42.50.Dv 



00 

o 
o 

o 

Q 



Oh; 

a ■ 

a ■ 

^ : 



> 

m 
p 

o 

oo 

o 



X 



In recent years impressive efforts have been made to 
quantify the entanglement of two-party states of quan- 
tum systems. This trend was highly stimulated by the in- 
terest in exploiting entanglement as an efficient resource 
in quantum information processing. For any pure bi- 
partite state a convenient measure of entanglement is 
now unanimously admitted, namely, the von Neumann 
entropy of its reduced states P, Unlike the pure- 
state case, several measures of entanglement have been 
considered for mixed bipartite states on both ffnite- and 
infinite-dimensional Hilbert spaces Q . Because of its op- 
erational meaning, the entanglement of formation (EF) 
of a mixed bipartite state, introduced by Bennett et al. 

plays a significant role: it is the minimal amount of 
entanglement of any ensemble of pure bipartite states re- 
alizing the given state. To be explicit, the EF of a mixed 
bipartite state p is defined as an infimum taken over all 
its pure-state convex decompositions 4]: 

Ep{p) :=inf{^pfc£;(|vl/,.)(vl/,|) |p = ^p,.|vl;,)(vl/,|}.(l) 

k fc 

Here E{Y^ k) k\) is the amount of entanglement of the 
pure bipartite state l^'fc). According to definition ([1]), 
evaluating the EF is a hard task, even for special quan- 
tum states. However, analytic evaluations of the EF have 
been carried out in a few finite-dimensional cases: gen- 
eral two-qubit states [H] , isotropic states [6] , and Werner 
states 0. 

In quantum information with continuous variables, 
two-mode Gaussian states (TMGSs) of the quantum ra- 
diation field are especially accessible from both theoret- 
ical and experimental standpoints. Their usefulness was 
recently reviewed in Refs. 8, 9]. So far, the only evalu- 
ation of the exact EF in an infinite-dimensional Hilbert 
space has been performed for symmetric TMGSs 
Moreover, the additivity of the EF has been proven in 
this case [ll|. The Peres-Simon separability theorem 
made it possible to use Gaussian measures of en- 
tanglement. Within a Gaussian approach, the reference 
set of states involved in the definition of any accepted 
entanglement measure is restricted to the subset of the 
Gaussian ones. Thus, following an earlier distance-type 
proposal for quantifying entanglement due to Vedral and 
co-workers [2|, several Gaussian evaluations employing 



the relative entropy 13 1 or the Bures metric [isl [l6l| 
have been performed. In Ref. a Gaussian entangle- 
ment of formation (GEF) has been introduced for any 
inseparable TMGS by analyzing its optimal decomposi- 
tion into pure TMGSs. 

The aim of the present work is threefold. First, we 
build the appropriate decomposition, Eq. ((T|), of a TMGS 
Pg that allows us to show that its EF and GEF coincide. 
We thus answer an op en p roblem in continuous-variable 
quantum information [lOL . Second, we give a more 
comprehensible approach to the problem of evaluating 
the GEF by use of covariance matrices (CMs). This en- 
ables us to write equations that yield, via the resulting 
optimal decomposition, an analytic solution for the EF in 
the general case. We also get explicit results in the most 
interesting special cases. Third, based on this approach, 
we prove the additivity of the EF for two-mode Gaussian 
states. 

Before proceeding we recall several useful properties 
of TMGSs. For later convenience, we choose to describe 
any TMGS pc by its characteristic function (CF), 



Xg(Ai,A2) := Tr[pG^i(Ai)I?2(A2)], 



(2) 



where D{a) := exp (aa^ — a* a) is a Weyl displacement 
operator. The CF of an undisplaced TMGS is xg{x) = 
exp (— ix-^Vx) . Here a: G R"* and V is the real, sym- 
metric, and positive 4x4 CM that completely describes 
the state. Its entries are the second-order moments of 
the canonical operators qj = {aj -\- aj)/-\/2, Pj = {oj — 

aj)/(-\/2«), where Uj and a^, (j = 1, 2), are the amplitude 
operators of the modes. Note that V e M4(R) is the CM 
of a TMGS if and only if the Robertson-Schrodinger ma- 
trix inequality holds: V -I- > 0, := i{a2 ® CT2), with 
0-2 a Paufi matrix. In particular, V :— det (V -|- > 0. 
Gaussian states whose CMs are connected by local sym- 
plectic transformations have the same amount of entan- 
glement and belong to an equivalence class: their CMs 
are locally congruent to CMs having a scaled standard 
form, 



V{ui,U2) 



( h\U\ Cy/uiU2 \ 

bi/ui d/y/uiU2 

C^MlU2 62W2 

y d/y/uiU2 &2/U2 / 



•(3) 
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In Eq. ([3]), ui > 1,U2 > 1 are one-mode squeezing factors. 
The unsealed standard form V(l,l) of the CM, intro- 
duced in Ref. , is expressed in terms of four parameters 
bi, 62, c, d. They are local invariants and determine the 
entanglement properties of the whole equivalence class. 
Recall the locally invariant Peres-Simon separability con- 
dition for a TMGS [H, V > 0, with V denoting the 
CM of the partially transposed density operator. This 
matrix inequality reduces to the Simon separability test 

V-.^det (j^ + ^nj ^ detV-^{bj + bl + 2c\d\) + ^ > 0. 

(4) 

The concept of classicality (existence of the Glauber- 
Sudarshan P representation of the density operator) is 
central in our present treatment of the EF. A TMGS 
with a CM ([3]) is classical if and only if the matrix 
V{ui,U2) — ^li is non-negative, with /4 the 4x4 iden- 
tity matrix. This requirement is equivalent to the non- 
negativity of all its principal minors. Remark that the 
classicality conditions are not locally invariant, depend- 
ing on the factors ui, U2. 

We start on the programme of Eq. ([1]) for an in- 
separable mixed TMGS pc, whose CM has a scaled 
standard form ([3]). Its four standard-form parameters 
&i,62,c > |d| — —d > are given, while the scaling 
factors ui , U2 are unknown. Continuous pure-state de- 
compositions of such a state are convex combinations of 
the type 

PG = yd2/3id2/32P(/3i,/32)|vI'(/3i,/32))(*(/3i,/32)|. (5) 

P{Pi, (32) is a non-negative normalized distribution func- 
tion and |\l/(/3i, /32)) is a state vector depending on the 
complex variables Pi, (32- In accordance with the EF def- 
inition, Eq. ID), the pure states in the above continuous 
combination should achieve an optimal decomposition of 
the given state pc- To this end, we make use of an im- 
portant theorem regarding the ranking of entanglement 
among pure two-mode states proven in a recent paper of 
Giedke et al. [l3| : For a given EPR uncertainty, the min- 
imal entanglement over the whole class of pure states is 
reached by a Gaussian one, the two-mode squeezed vac- 
uum state (TMSVS). This important result leads to the 
key idea of our treatment: Owing to the Gaussian nature 
of the two- mode state pa, as well as to the scaled stan- 
dard form ^ of its CM, we are allowed from the very 
beginning to restrict ourselves in Eq. (O to equally entan- 
gled pure states obtained by displacing a unique TMSVS. 
Among all ensembles of such pure two-mode states that 
realize the given mixed state pc we have to find the one 
possessing the minimal entanglement. Let us denote by 
Pq = the TMSVS entering this optimal convex 



expansion: 

PG^ J d^pid^P2Pif3i,P2)Di(f3i)D2iP2)poDl{l32)Dl{/3i). 

(6) 

According to Eq. ([6]), the exact EF of the given mixed 
two-mode state pc reduces to the amount of entangle- 
ment of the TMSVS po: 

Ef{pg) = E{po). (7) 

Recall now that a TMSVS is a Gaussian state whose CM 
has precisely the unsealed standard form V(l, 1), Eq. ([3]): 
its parameters bi — b2 —'■ x > 1/2, c = —d =: y > are 
subjected to the purity condition 

^'-y' = l- (8) 

The entanglement of a TMSVS is the von Neumann en- 
tropy of its one-mode reductions, 

EiPo) = i^+l) Mx + l)~ix-l) ln(x - i), (9) 

which is an increasing and concave function of the vari- 
able X > ^. We notice that the optimal convex de- 
composition ([6]) stands the pertinent test of the pure- 
state limit case: pa = Po, with ui = U2 = 1 and 
P{Pi,l32) = S2{Pi)S2{f32). To evaluate the EF jT]), one 
has to be able to determine the optimal decomposition 
i.e., both the distribution hmction P(/3i, ^2) and the 
TMSVS Po- Provided that this can be effectively done 
for any TMGS, Eq. ([6]) displays the first result of our 
work: the EF for a two-mode Gaussian state coincides 
with its GEE. 

We now take advantage of a fact well known in 
quantum optics: decompositions of the type ^ do 
have a clear meaning starting with Glauber's semi- 
nal work on the coherent states of the electromag- 
netic field pJi]. Accordingly, Eq. ([6]) gives the den- 
sity operator pc of a superposition of two fields: one 
is in a classical state pc having the regular Glauber- 
Sudarshan P representation P(/3i,/32), and the other is 
in the pure state po. By employing the corresponding 
CFs and writing the P representation as the Fourier 
transform of the normally-ordered CF, x'^^^(Ai,A2) := 
exp(i(|Aip + |A2p))x(Ai, A2), Eq. ® leads to the mul- 
tiplication law 

x[f^(Ai,A2) =xr(Ai,A2)x[f^(Ai,A2), (10) 

with xc denoting the CF of the classical state pc of 
the superposed field [3, [2^. It follows that the CF 
Xc(Ai,A2) is also Gaussian. Equation (fTO| results in 
an addition rule for the CMs of the Gaussian states in- 
volved: 

V{ui,U2) = Vo + Vc-\h. (11) 
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Our method towards finding the optimal pure-state de- 
composition concentrates on the properties of the classi- 
cal state PC- We first show that pc belongs to the bound- 
ary dV of the set V of all classical TMGSs [Property 1). 
Then we prove that pc is also on the boundary dS of the 
larger set S of all separable TMGSs: S D V [Property 



A2 



2), Refs. [21|,|22| 



Property 1: The superposed classical state pc is at the 
classicality threshold. Indeed, for the optimal superpo- 
sition, the CM V[ui,u2) should be as close as possible 
to Vq. This happens when the principal minors of rank 
3 and 4 of the non-negative matrix V(ui,U2) — Vq = 
Vc — \li are zero. By the same token, the Gaussian 
state Pc is at the border of classicality &P . Explicitly, 
the condition det(Vc — ^14) = holds with the left-hand 
side expressed as a product of two vanishing factors: 

- x)[b2U2 ~x)- [cy/uiU2 - yf = 0, (12) 



- x){h2/u2 - x) - (Ml/V^ - V? - 0. (13) 



Equations (|T2|) and (|T3|) are in agreement with the 
Gaussian optimality conditions written in the pioneering 
work Ref. ll| on different grounds. Making use of Eqs. 
([HI), p2)l . and ((T3| . we can impose to the one- variable 
function x = x[ui,U2[ui)) the minimization condition 
= 0. We get therefore a fourth independent algebraic 
equation. 



biux 



b2lu2 



(14) 



which implies an additional property of the state pc- 

Property 2: The superposed classical state pc is at the 
separability limit as well. To prove this statement, we 
use Eqs. ([5]), (HI]), and (fT5)) to evaluate the Simon in- 
variant V, Eq. (dl, of the Gaussian state pc- Taking 
into account Eq. (fTi)) . we get V = 0, i. e., pc G dS. 

The evaluation of the required EE reduces to solving 
a system of four non-linear algebraic equations, namely, 
Eqs. and (fT^ - (fH|) . with four unknowns: Mi, 1*2, x, y. 
Let us denote its solution by wi,W2,Xm,ym- The above 
algebraic system yields a quartic equation, J2n=o AnP" = 
0, for the product p :— uiU2- The coefficients An are 
quite simple polynomials in the four standard-form pa- 
rameters of the given inseparable TMGS: 



Ao = (6162 -d') 
X b2[bib2 - 



biibib2 - d^) - ^ 



>0, 



Ai = -[c(6i62"d') + ^]|(6i-&2)'[c(6i62-d') + ^] 



-26162(0- 



6162 - d^ - i 1 < 



[(6ic-62|d|)(6i|d| 
-2[blbl - c^d^)V - 



-62c)-hc|d|2:](detV-l-l/16) 
c|d|detV, 



^3 ^^i(c^ 



A(c ^ Ml) > 0. (15) 



We have introduced the symplectic invariant Z :~ b\-\- 
62 + 2cd >l/2. Atp=l, the above quartic polynomial 
has a negative value, except for c = |c?|, when it vanishes. 
This implies the existence of a convenient root pm = 
W1W2 > 1 for any inseparable mixed TMGS. Had we got 
Pm, it could be used to obtain the optimal as the 
smallest root of a quadratic trinomial B2[p)y^ + Bi[p)y + 
Bo{p) whose coefficients, 

Boip) = -Vp > 0, 

Bi[p) = -2^[[\d\[bib2-c^) + c/4]p 
+ [c(6i62-d') + M|/4]) <0, 

B2[p) = (6162 - c2)p2 + Zp+ (6162 - d^) > 0, (16) 

are evaluated at p — pm. We mention that in four sig- 
nificant particular cases (defined by special relations be- 
tween standard-form parameters) we have found simple 
solutions by direct use of Eqs. (fT2|) - ((T4)) . We have then 
recovered them by exploiting Eqs. (fT3|) and (fTB)) . 

As a first salient example, we consider an entangled 
symmetric TMGS, whose standard- form parameters are 
61 — b2 —■ b, c> \d\ = —d > 0. The smallest symplec- 
tic eigenvalue k_ of the CM for the partially transposed 
density operator is in this case k_ = ■\/(6 — c)(6 — |d|). 
In agreement with the results of the remarkable work 
Ref.[l3], Eqs. ^ and §^ give: 



Wi ^ W2 = 



b-\d\ 



1/4 



2k_ 



A second class of notable bipartite states is that of 
two-mode squeezed thermal states. The standard-form 
parameters of such a state are 61 > 62, c — \d\ = —d > 0. 



This case was considered previously in Refs. [23|, |2J], 
where the prescription of Ref.[ll[ to evaluate the GEF 
was followed. From our results, 



Wi = W2 = 1, x„ 



(61 + 62)(6i62 - 1/4) - 2cVV 
(61 + 62)2 - 4c2 

is not determined only by the eigen- 



one can see that a;„ 
value K_ = i[6i -f 62 - V(6i - 62)^ +40^]. 

A third example is that of a TMGS at the separability 
boundary: V — <J=4> — ^. We get Xm = ^, ym = 
and the optimal squeeze factors 



Wi = 



Ib2[bib2-d^)-\bi 

b2{bib2-c^)-lbi' 



W2 — Wi(6i ^ 62). 
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Equation ([6]) becomes now the P representation of a 
state pg at the border of classicahty and that of sep- 
arabihty dS as well. 

A fourth class of entangled states consists of those 
TMGSs whose CMs have the smallest symplectic eigen- 
value K_: I? — ■<F=^ K_ = ^. These states were studied 
as having minimal negativity at fixed local and global 
purities [25|]. Assuming that hi > 62, c > \d\ — —d > 0, 
we found two distinct solutions required by the sign of 
the difference 620 — 5i|c?|: 



Wi 



62C — 6i|(i| < : x„ 



Ib2{blb2-d^)-\h 
62(fol&2-c2)-i6i 



62C — bi\d\ > : x„ 



bl 



8(detV-^) 



W2 = Wl(&l ^ 62)- 



6162 



6i&2-rf2' 



Wi 



bl 



(6162 — d^), 'W2 — 'Wi{bi ^ b2) ■ 



The above formulae for Xm are in agreement with those 
derived in other parametrization in Ref. f24!|, which fol- 
lows the methods of Ref. 11 1 . 

The last issue we are here interested in is the additiv- 
ity of the EF for TMGSs. Our present approach gives 
a straightforward answer to this open question We 
consider a four- mode product state pc ® cgj where pc 
and aa are entangled TMGSs. We denote the minimally 
entangled TMSVSs entering the optimal decompositions 
of the type ^ for both factors by po and ctq , respectively. 
Therefore, their tensor product pqiExtq enters the optimal 
convex decomposition of the four- mode state pg®(^g- It 
follows the identity EpipG^crG) = E{poiS)ao). The well- 
known additivity property of the von Neumann entropy, 
E{po (8) (To) = E{po) + E{ao), yields the additivity of the 
EF for TMGSs: 

Ef{pg «> ctg) = Ef{pg) + Ef{(Jg)- 

Consequences of this property on evaluating other mea- 
sures of entanglement are largely discussed in Ref. 3] . 

To sum up, we have reformulated the problem of eval- 
uating the EF for TMGSs in terms of CFs and CMs. We 
have shown that the exact EF of such a state coincides 
with its Gaussian one. Although an analytic solution in 
the general case seems to be complicated, it can be found, 
nevertheless, by solving a quartic equation. Our general 
treatment allowed us to retrieve readily previous explicit 
results in some relevant particular cases. Based on our 
approach, we have finally proven the additivity of the EF 
for two-mode Gaussian states. 
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Note adrferf.-During the completion of this Letter, an 
interesting treatment of the EF for a TMGS was given in 
Ref. [2^ . Its relation to our present work will be discussed 
elsewhere. 
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